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Asymptotic properties of the solution of two-dimensional 
randomly forced Navier-Stokes equation with long-range cor- 
relations of the driving force are analyzed in the two-loop 
order of perturbation theory with the use of renormalization 
group. Kolmogorov constant of the energy spectrum is cal- 
culated for both the inverse energy cascade and the direct 
enstrophy cascade in the second order of the e expansion. 

47.27.Gs,05.10.Cc 



I. INTRODUCTION 

Asymptotic properties of the stochastic Navier-Stokes 
equation have been actively studied with the aid of renor- 
malization group (RG) during the last two decades (for 
a comprehensive review, see A random force with 
powerlike falloff of spatial correlations has been used to 
maintain a stationary scale-invariant regime with a sub- 
sequent generation of an e expansion of scaling exponents 
and scaling functions. 

Most work has been carried out in view of appli- 
cation to three-dimensional turbulence; partly because 
of greater physical interest than in the two-dimensional 
case, partly due to additional divergences, which occur 
in two dimensions and prevent a direct use of the previ- 
ously obtained general d-dimensional results. Moreover, 
there were some serious flaws in the early work |^ on 
the RG-analysis of the forced Navier-Stokes equation in 
two dimensions, and only recently has a consistent renor- 
malization procedure been put forward |^,| | and used in 
various problems involving randomly forced incompress- 
ible fluid 

Apart from energy, in two dimensions the enstro- 
phy (squared vorticity) is an inviscid conserved quan- 
tity quadratic in the velocity. Therefore, two self-similar 
regimes corresponding to an inverse energy cascade to- 
wards small wave numbers and a direct enstrophy cas- 
cade towards large wave numbers are expected to take 
place 1^ instead of the direct energy cascade observed 
in three-dimensional turbulence. The energy (enstro- 
phy) pumping leading to a steady state with the two 
scaling regimes may be realized in two different ways. 
On one hand, in numerical simulations ||^ and some ex- 
periments the energy and enstrophy pumping takes 
place on scales in between the inverse energy cascade 
and the enstrophy cascade. On the other hand, in atmo- 



spheric turbulence pl] | the energy and enstrophy sources 
are at the outer edges of the scaling intervals, and it 
is not clear whether there is an energy and enstrophy 
sink between them or they coexist Q . In both cases 
the Kolmogorov spectrum of the inverse energy cascade 
E{k) cx k^^/'^ for fc <C fc/ is observed experimentally and 
in the majority of simulations for wave numbers smaller 
than the characteristic wave number kj of the energy 
pumping. For the enstrophy inertial range the existing 
data is not so clear-cut, although recent high-precision 
numerical simulations [|2| and accurate analysis of the 
atmospheric data seem to support the existence of 
an energy spectrum cx k^^ in the enstrophy cascade. 

In this paper we have carried out a two-loop renor- 
malization of the field theory generated by the two- 
dimensional randomly forced Navier-Stokes equation. 
We have taken into account the divergences specific of 
two dimensions in the way proposed in Ref. |^], and cal- 
culated the asymptotic expression for the energy spec- 
trum in the second order of the e expansion for both the 
energy and enstrophy inertial ranges. For convenience of 
calculations, we have used the stochastic Navier-Stokes 
equation with a subsequent dimensional regularization of 
divergences in two dimensions Q instead of the specifi- 
cally two-dimensional setup Q with the stream-function 
description of the fiuid flow. 



II. RENORMALIZED FIELD THEORY FOR THE 
TWO-DIMENSIONAL STOCHASTIC 
NAVIER-STOKES EQUATION 

Consider the stochastic Navier-Stokes equation for the 
flow of homogeneous incompressible fluid, which for the 
transverse components of the velocity field assumes the 
form 



dtu, + PijVidiVj = t'oV^Wi - ^oVi + Ft 



(2.1) 



together with the incompressibility condition diVi — 0. 
In Eq. (2T) Vi{t,x.) are the coordinates of the diver- 
genceless velocity field, vq is the kinematic viscosity, ^0 
is the coefficient of friction, and Py is the transverse pro- 
jection operator (P^ = 6ij — kikj/k'^ in the wave-number 
space), and Fi are the coordinates of the random force. 
Here, and henceforth, summation over repeated indices 
is implied. 

In experimental realizations and simulations of a two- 
dimensional turbulent flow energy may be consumed not 
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only by microscale dissipation, but also by the friction 
at th e bo undaries of the fluid layer. The friction term in 
Eq. (2.1) makes it possible to maintain stationary state 
with the anticipated inverse energy cascade towards small 
wave numbers and the direct enstrophy cascade towards 
large wave numbers, when the pumping is carried out in 
between the corresponding inertial ranges. The coeffi- 
cient of friction is a mass term from the point of view of 
renormalization, therefore we put = in the calcula- 
tion of the renormalization constants of the solution of 
Eq. (pj|). As shown in Rcf. [Q, the friction term in Eq. 
( ^.l[ ) is not renormalized and thus does not affect the RG 
equations and the subsequent asymptotic analysis. 

In the applications of the stochastic Navier-Stokes 
equation (2T) to turbulence the random force is assumed 
to have a gaussian distribution with zero mean and the 
correlation function in the wave-vector space of the 
form 

(F, (i, k)F, (f, k')) = P^J {2nf^''5{k + k')6{t~t')dF (fc) . 

(2.2) 

The scalar kernel has a powerlike asymptotic behavior at 
large wave numbers: 



dpik) = Dok 



4-d-2e 



h{m/k) 



(2.3) 



where h{x) is a well-behaved function of the dimension- 
less argument m/k ensuring the convergence of the in- 
verse Fourier transform of dp{k) at small k and with the 
large k behavior fixed by the condition h{0) — 1. In a 
fixed dimension above two dimensions the force correla- 
tion function is not renormalized and the the kernel ( |2.3| ) 
remains intact. This is, however, not the case in two di- 
mensions, in which renormalization generates additional 
terms oc k^ into the force correlation function. In order 
to deal with a multiplicatively renormalizable theory - 
which is convenient technically - we add this term to the 
correlation kernel at the outset and use, instead of the 



function (2.3), the modified function 



d'p{k) = Doik 



h{m/k)+Do2k^ 



(2.4) 



The force correlation function is related to two basic 
physical quantities, the energy pumping rate £ and the 
enstrophy pumping rate B, as 



£ = 



B 



d-1 



dk 



{2n) 
dk 



(2^) 



(2.5) 



in d-dimcnsional space, which allows to connect the " cou- 
pling constant" Z?oi with the pumping rate in the corre- 
sponding asymptotic region. 



We cast the stochastic problem (2.1), (2.2), (2.4) into 



a field theory with the De-Dominicis- Janssen " action" in 
the usual manner . An analysis of UV divergences with 



the use of Galilei invariance, causality and symmetries of 
the model allows to write the renormalized action in the 
form 



1 



dtdk 
{2'kY 



giv^l?%k^Y-^-^h{m/k) + g2iy^n^^^Z2k^ 



(2.6) 



+ J dtdxv ■ [dtv + (v • V)v - i/ZiV^v] , 

where is the scale-setting parameter of the renormal- 
ized model and 25 = d~2 is the parameter of dimensional 
regularization. Only two renormalization constants Zi 
and Z2 are needed to absorb the UV divergences of the 
model in two dimensions. To avoid excessive notation, 
we have used the same symbols for both the fields and 
their Fourier transforms in (p.6|). 



Renormalized parameters of the action (2.6) are de- 
fined by 



i^o = i^Zi , 
D02 = .g2i^V"'''^2^r' 



(2.7) 
(2.8) 
(2.9) 



We have used a combination of dimensional and analytic 
regularization with the parameters e and 2(5 = rf — 2. As 
a consequence, the UV divergences appear as poles in 
linear combinations of the regularizing parameters. Nor- 
malization has been fixed by the choice of the minimal 
subtraction (MS) scheme jl^ . 



III. RENORMALIZATION-GROUP EQUATIONS 
AND FIXED POINTS 

We set up the notation and basic equations for the 
spatial Fourier transform of the pair correlation function 
of the random velocity field 



Wrnnih - t2,k;Doi,DQ2,iyo) = J 



dV 

{2ttY 

x(^;™(^l,Xl)«„(^2,X2))e^''•(^l-^^^ 



(3.1) 



because this quantity is directly connected to the en- 
ergy spectrum through the relation (u„(<, x)D„(t, x)) = 
2j^°°E{k)dk. 

Independence of the unrenormalized pair correlation 
function 

Wrnn{t,k; Doi, 002,1^0) = W^„{t,k; gi, g2, 1^, fJ.) 

of the velocity field v of the scale-setting parameter /i 
gives rise to the basic RG equation 



9 n 9 ^ d d 
^^l dgi dg2 ov 



= (3.2) 
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for the renormalized cor rela tion function Wj^„ . The coef- 
ficient functions of Eq. ( |3.2| ) Pi, (32, and 71 are expressed 
in terms of logarithmic derivatives of the renormalization 
constants. We use the definitions 



dlnZi 



(3.3) 



where i = 1,2, and the subscript "0" refers to partial 
derivatives taken at fixed values of the bare parameters. 

It is convenient to express the correlation function 
through a dimensionless scalar function R. In two- 
dimensional space we define this function through the 
relation 

W'm«(i,k;gi,52,J^) = i5ii'^-Pmn(k)i?(T, s; gi, 32) , (3.4) 

where s — fc//i is the dimensionless wave number, an d 
T = ivkl^ the dimensionless time. Solving Eqs. ( ^.2[ ), 
( |3.3|) by the method of characteristics we obtain the cor- 
relation function in the form 

W^„(t,k;5i,g2,i') = ^giV^P^n{k)R{tk^V,l;gi,g^) , 

(3.5) 

where are the solution of the Gell-Mann-Low equa- 
tions: 



din s 



Pi [91,92] 



1,2. 



and v is the running coefRcient of viscosity 



(3.6) 



(3.7) 



Writing the latter in terms of the unrenormalized (physi- 
cal) parameters and the running coupling constant 'g^ |^] 
as 



Dpi 

9i 

we arrive at the expression 



1/3 



(3.8) 



W^^{t,k;gi,g2,iy) = -f/'Dli'k-'^^' 
xP„„(k)i? ((i?oi/5i)'/'fc'-''/'t, l;5i,52) ■ 
For the /3-functions: 

/3i = .gi(-2£ + 37i), 
P2 = .92(^72 + 371) , 



(3.9) 



(3.10) 



a tedious two-loop calculation, with the use o f th e step 
function h{m/k) = 9{k — m) in the kernel (2.4), and 
dimensional rcgularization, yields 



71 = iLi +U2 + Ui y-a + 



72 



[ui + U2Y , w?(l - r) 



U2 



U2 



(3.11) 
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,'5 25 \ 2/1 
+W1M2 [ ^a + — -irj +U2 I -- - r 



where 



Ul 



91 
327r 



U2 



92 
327r 



(3.12) 



The constant r = —0.1685 in Eq. ( 3.11 ) comes from a 
numerical calculation of those parts of two-loop graphs, 
for which analytic results were not feasible. The con- 
stant a = C - ln(47r) = -1.9538 (C is Euler's con- 
stant) is brought about by a 2i5 = c? — 2 expansion of 
the geometric factor Sd/i2nY = 2/[r(d/2)(47r)'*/2] = 
^[1 + aS + 0{6^)]. The method of calculation is es- 
sentially the same as was used for the d-dimensional case 
in Ref. ITll 



At one- loop order the 7- functions (3.11) are exactly the 



same as those of the d-dimensional Navier-Stokes equa- 
tion in two dimensions |3|. They also coincide with the 
expressions obtained directly in two dimensions for the 
corresponding stochastic vorticity equation Thus, we 
think that for calculation of the coefficient functions of 
the renormalization group equation the results of the d- 
dimensional model in the two-parameter expansion may 
be applied directly to the two-dimensional case. There 
might be some discrepancies in calculations involving 
composite operators due to different symmetries in two- 
dimensional and general d-dimensional cases, but we do 
no calculate anything like that here. 

The fixed points are determined by the syste m of equa- 
tions Pi = P2 = 0. From the solution of Eqs. (3.6) near 
a fixed point it follows that the fixed point is infrared 
stable, when the matrix — dnPm is positive defi- 
nite. If £ < 0, then the trivial fixed point: = U2 = 
is infrared stable. The anomalous asymptotic behavior 
of the model at small wave numbers is governed by the 
nontrivial fixed point 



4 = le_ A(2a + 5-4r), 



9 27 
2 4 

" 9^~ 81 



(a - 2 - 3r) , 



(3.13) 



at which the eigenvalues of the stability matrix are 




(3.14) 



3 



Note that these eigenvalues a re ind ependent of a. The 
real parts of both eigenvalues ( 3.14 ) are positive and the 
fixed point (3.13) infrared stable, when e > 0. 



IV. TWO-LOOP CALCULATION OF 
KOLMOGOROV CONSTANTS 

The connection between the energy spectrum -E(fc) and 
the equal-time correlation function of the velocity field, 
(w„(x)w„(x)) = 2 E{k)dk, in the two-dimensional 
wave- vector space amounts to 



E{k) = AT4.„„(o,k) 
47r 



From (3.9) the asymptotic expression 



jLl-4£/3 



(4.1) 



(4.2) 



follows, when fc ^ 0. Here, gl — 32ttuI, — 32nu2 are 
the values o f the coupling constants at the infrared-stable 
fixed point ( 3.13 ). 

The relations (2.5) allow to express the parameters Dqi 
and Dq2 in terms of the energy (or enstrophy) pumping 
rates £ (B) . Integrating over the wave- number shell m < 
k < A with h{x) = 1 in the kernel ( ^.4| ) we obtain, in the 
limit of widely separated upper and lower wave-number 
limits. 



Doi A2(2-^) 



B 



8n 2-e 

Dpi A2(3-e) 

"8^ 3-e 



02 

IGtt 
247r' 



A'- 



(4.3) 
(4.4) 



The spectrum ( p~^ ) should be independent of the details 
of the energy pumping, i.e. independent of the upper cut- 
off A in the range m <C fc ^ A. According to the relation 
(4.3), this goal is achieved by the choice D02 = and 
e = 2 for the anticipated inverse energy cascade. The 
relation (4.4), in turn, shows that the choice D02 = 
and £ = 3 leads to scale-invariant behavior for the di- 
rect enstrophy cascade. In both cases it should be borne 
in mind that the bare coupling constant Dq2 is techni- 
cally a book-keeping parameter refiecting the necessity 
of the introduction of the short-range term in the corre- 
lation function of the random force. Physically, it could 
be related to the intensity of thermal fluctuations, which, 
however, are irrelevant in the energy balance of station- 
ary developed turbulence. 

The Kolmogorov constants are determined from the 
asymptotic relations 



E{k) = C(£)£2/3fc-5/3 (^^^ 



E{k) = C'(e)S2/3fc-3 



A 



4(e-2)/3 



4(£-3)/3 



(4.5) 
(4.6) 



Thus, the choice e = 2 in (4.3) renders the spectrum ( |4.2| ) 
completely scale-invariant with the Kolmogorov expo- 
nents corresponding to the energy cascade, whereas the 
substitution e = 3 in (4.4) leads to scale-invariant behav- 
ior in the enstrophy cascade. 

The Kolmogorov constants may be calculated in the e 
expansion from (4.2), (4.5) and ( |4.6| ). The present two- 
loop calculation allows to find correction terms to the 
previosly found B expressions. The result is 



C{e) = 2 ■ 3i/3£i/3 
C'{e) = 3 • 2^^^e'/' 



l + -{l + r)e 



' 3 + 2r ' 
1 + — z — e 
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(4.7) 
(4.8) 



For e = 2 we obtain from (^ C = 4.977. The clo- 
sure model leads to the prediction C ~ 6.69 Re- 
sults of numerical simulations vary from C — 2.9 |l6|| to 
C 14 [|l^. Experimental results jlO| yield the range 
3 < C < 7. The result obtained here is thus in better 
agreement with other available data than the leading- 
order value of the e expansion C = 3.634 obtained in 
Ref. |. 

For e = 3 the value C = 10.29 obtained from (HJ) 
is significantly larger than the leading-order value C — 
5.451 of Ref. [|j and the discrepancy between the present 
result and the closure-model prediction C" = 2.626 ||] 
is larger. However, from a more detailed analysis of the 
model it may be concluded that calculation of the con- 
stant C is not unambigous. The value obtained from 



(4.8) corresponds to the case, in which the coefficient of 
friction ^0 = 0- If? however, this coefficient is retained, 
then the asymptotic expression for the energy spectrum 



xR 0,l;gl,g^, 



1/3 

9* i^m 

\l/3 



Mr!eo,-(2- 



1/3 
01 



87r 

2e/3) ] 



(4.9) 



instead of (4.2). The energy spectrum ( [O] ) is scale- 
invariant for e = 3 regardless of the value ^o- Therefore, 
the proportionality constant in the scaling law seems to 
be nonuniversal in the enstrophy inertial range, and de- 
pends on the properties of large-scale dissipation Q . Re- 
cent spectral closure analysis and numerical simulations 
have led to similar conclusions |ll8| . 



V. CONCLUSION 

In this paper we have carried out a two-loop renormal- 
ization of the randomly forced Navier-Stokes equation 
with long-range correlated random force in two dimen- 
sions in view of two different patterns of scale-invariant 
asymptotic behavior. 
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We have calculated the Kolmogorov constant for a 
powerlike asymptotic energy spectrum cx k~^^^ of the 
random velocity field in the inertial range of the inverse 
energy cascade in the second order of an e expansion with 
the result C = 4.977, which is in reasonable agreement 
with other available experimental and theoretical data. 

We have also calculated the Kolmogorov constant for 
the spectrum oc in the inertial range of the direct en- 
strophy cascade with the second-order result C = 10.29 
with a larger deviation than at the leading order from 
results obtained by other methods. However, explicit 
asymptotic expressions for the pair correlation function 
of the random velocity field obtained in the present ap- 
proach strongly indicate that the Kolmogorov constant 
in the enstrophy cascade is not universal, but depends on 
the enstrophy dissipation due to large-scale friction. 
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